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the Deutsch problem 

M. M. Santos, F. 0. Prado, H. S. Borges, A. M. Alcalde, J. M. Villas-Boas, and E. I. Duzzioni 
Universidade Federal de Uberlandia, Caixa Postal 593, 38400-902, Uberlandia, MG, Brazil 

Abstract 

The wide set of control parameters and reduced size scale make semiconductor quantum dots attractive candidates to 
implement solid-state quantum computation. Considering an asymmetric double quantum dot coupled by tunneling, we combine 
the action of a laser field and the spontaneous emission of the excitonic state to protect an arbitrary superposition state of the 
indirect exciton and ground state. As a by-product we show how to use the protected state to solve the Deutsch problem. 



I. INTRODUCTION 

One important feature of Quantum Mechanics is the superposition principle. This fundamental characteristic for 
pure states can be summarized by the quantum entanglement phenomenon [l], Q • It is believed that such a property 
can give to the quantum world the advantage of processing information in a more efficient way than its classical 
counterpart Q. However, in a realistic scenario, decoherence destroys the superposition of states leading the system 
state to a statistical mixture |4|. A number of methods has been proposed to circumvent this difficulty; among 
them we highlight the quantum error-correcting codes [5|, decoherence- free subspaces in collective systems |6j, and 
dynamical decoupling methods [7j. Differently from the methods cited above, the engineering reservoir technique 
[8| makes use of incoherent control of Markovian reservoirs to drive the system state to a desired superposition of 
pure states in the asymptotic limit [9]. Therefore, the effective interaction between the system and reservoir must 
be carefully engineered to drive the system to equilibrium with the reservoir. Recently this theme received a great 
deal of attention in order to obtain pure entangled protected states Q]. Besides the recent generalization of the 
engineering reservoir technique to protect nonstationary superposition states [lOMl2l . two important applications of 
this technique have emerged: (i) the construction of robust quantum memories }13| and (ii) the implementation of 
quantum computation via dissipation (QCD) 14]. In the former case the quantum information can be stored for long 
times at the nodes of a quantum network, while in the later case the quantum processor is insensitive to external 
perturbations or decoherence. Particularly, in Ref. [14] it was shown that QCD is universal and can efficiently simulate 
a given quantum circuit. 

One of the clearest demonstration of the power of quantum processing is the quantum solution for the Deutsch 
problem [15J. Later improved by Deutsch and Jozsa [16J ] and by Collins, Kim, and Holton [17| . the optimized 
version of the Deutsch algorithm is able to decide wheter a binary function is constant or balanced with a single 
measurement, while in the classical case two measurements are required. The Deutsch problem was proved using 
different experimental setups, such as Nuclear Magnetic Resonance [T|| Ha], optical systems [2Q], circuit quantum 
electrodynamics architectures [2l| , trapped ions [22[ , the nitrogen- vacancy defect center [23| , and quantum dots [24| . 

In this work we investigate the protection of superposition states in a physical system composed by two asymmetric 
quantum dots coupled by tunneling under the application of laser fields, where the interplay of the optical field and the 
tunneling creates a structure of two excitonic levels: direct and indirect, as described in Ref. [25| . The results obtained 
here generalize the result obtained in Ref. [26|, where a robust state of the indirect exciton was found. Through the 
protection of an arbitrary coherent superposition of the ground and excited (indirect exciton) states, we propose 
the implementation of the optimized version of the Deutsch algorithm [17| using QCD. In general, the proposals 
for implementing quantum computing protocols make use of coherent control of quantum systems. In contrast, our 
proposal to solve the Deutsch problem has the advantage of being controlled incoherently by the reservoir, becoming 
immune to decoherence processes. 

The manuscript is divided as follows. In Sec. II we obtain analytically the protected asymptotic state of the system 
which is parametrized on the Bloch sphere. Such a state is used to solve the Deutsch problem in Sec. III. Conclusions 
and perspectives for this work are presented in Sec. IV. In the Appendix we deduce the effective master equation 
leading to the protected state. 

II. PROTECTED PURE STATE 

We consider two asymmetric quantum dots coupled by tunneling, which under the application of a static electric 
field along the growth direction can be simplified to a three- level system [25|]: the ground state |0) (no excitation on 
the system), the exciton state |1) (one electron-hole pair in the same dot), and the indirect exciton state |2) (the 
electron in one dot and the hole in another). The exciton state can be created from the ground state by the application 
of a laser field with right the frequency and the indirect exciton from the exciton state by tunneling of one electron, 
as experimentally demonstrated in Ref. [27] . In this way, the Hamiltonian in the rotating wave approximation for 
this system is (h = 1) 

2 

H (t) = J2"i Ii) 01 + T e (|1) (2| + |2) (1|) + n [e^+^ |0) (1| + e -^+^0 |1) (0|] , (1) 

where u)j is the energy of the jth level and T e is the tunneling coupling between the levels |1) and |2). We tune the 
gate voltage to have u>i = u)2 = ojl and luq = 0. ui and tp are the frequency and phase of the laser field, respectively. 
fl = (0| ~ft ■ E |1) is the dipole coupling between the excitonic transition, with ~jt being the electric dipole moment and 
E the incident electric field of the laser [28[ ■ According to Ref. [26[ , the dominant decoherence processes can effectively 



be described by spontaneous emission from states |1) to |0) and |2) to |0) with decay rates Ti and r 2 , respectively. 
This approximation remains valid if we consider that (i) the system is at very low temperatures so that the pure 
dephasing induced by acoustic phonons does not modify significantly the system dynamics |2J|; (ii) the creation of 
the lower energy exciton state is made by the application of a resonant laser field, which inhibits transitions mediated 
by phonons [30, [3l| ; (iii) the optical phonon effects are negligible because they have quite different frequencies (~ 30 
meV) compared with that of our system (see the main text) |29|]. Then, the dynamics of the system can be described 
by the master equation 

Ml = -i[H(t), P } + /:(p), (2) 

with H (t) being the Hamiltonian defined in Eq. ([T]) and C (p) the dissipative Liouvillian, given by 

2 

2 



C (p) = & (p) + £ 2 (^Ef ( 2 1°) (*l P l»> <°l - 9 \i) (i\ ~ \i) (i\ p) ■ (3) 



The steady states of the system can be found through the condition 

hm d Jp- = 0, (4) 

where t — > oo means i>l/ minTi. The right-hand side of Eq. © fulfills the condition (JU) for a pure state |\I/) (also 
called dark state) if C (|\I>) (*|) = and H |<F) = E |*) [9]. For a master equation in the Lindblad form pS^ . 

C(p) = 1 -(20pO^ - O^Op- pO^O) , 

where 7 is the decay rate and O the jump operator, the state |'F) is the only protected state if O |W) = (C (|^) (^>\) = 
0) and there is no further eigenstate \<j>) of O such that [O, Q T 1 \<j>) = [33|. The method used here to obtain the 
dark state of the system is not unique. For instance, in Refs. [3J, |35[ the dark state conditions are met when the 
absorption spectrum is null. Already in Ref. 36] a signature of the dark state is found in the second-order correlation 
function. In order to remove the time dependence of the right-hand side of Eq. ([2]), we move to a rotating frame 
defined by the unitary transformation 



U (t) = exp 



^(|l)(f| + |2)(2|-|0)(0|) 



(5) 



In this frame the Hamiltonian becomes 

H mt = n (e* |0) (1| + e- 1 * |1) (0|) + T e (|1) (2| + |2) (1|) , (6) 

the form of C (p) remains unchanged, and p is replaced by pi nt = U* pU . Considering the requirements to obtain a 
protected pure state, we initially find the eigenvectors of H m t 

|i. + (^))^ Sin(g/2)|0)+e ^ |1 / l +e ^ COS(g/2)|2) , (7a) 

v2 

\Eo (8, if)) = cos (9/2) |0) - e~ iv sin (6/2) |2) , (7b) 

\E- (0, V) ) = Bin(g/2)|0)- e -^|l) +e -^cos(g/2)|2) i (?c) 

v2 



with eigenvalues E± = ± \/fl 2 + T e 2 and Eq — 0. tp is the laser phase defined above and cos (0/2) — T e / \/£l 2 + Tf. 
Through the condition £(|^) (^1) = we observe that the dissipative Liouvillians C\ and £2 have {|0) , |2)} and 
{|0) , |1)} as their dark stat es, resp ectively. As the relation between the decay rates is about T 2 = I0~ 4 ri with Ti of 
the order of 0.33 — 6.6/zeV |37H39| . we conclude that the dissipative dynamics is basically governed by C\ (p). Note 
that the eigenvector \Eq (9,<p)) is composed only by {|0) , |2)} states. Therefore, \E (9,<p)) is a dark state of the 
system, which is in good agreement with our numerical calculations performed in the regime of parameters defined 
above. In the Appendix we analytically show how to obtain the effective master equation whose protected state is 



\Eq (9, tp)). Naturally the protected state is not pure; the deviation from \Eq (9, ip)) introduced by the indirect exciton 
decay channel (T2) can be obtained through the fidelity J- (00) in the stationary regime, defined by 

1 1 r 2 r 2 ^ 
J" (00)= lim [(E Q (9,p)\p(t)\E a (9,p))} = r Fl ^S ■ (8) 



f 



Ti T 2 (T 2 +tt 2 ) 



A simple analysis of the particular case T e ~ 57 shows that the fidelity J 7 (00) ~ 1 — I^/Fi attains values next to 1, 
even for T2 one order of magnitude lower than T\. 

The state \Eq (9, tp)) can be represented on the Bloch sphere, where 9 and p are polar and azimuthal angles. The 
experimentally accessible values of fl ~ 0.05 — f .0 meV [38|, |4(| and T e ~ 0.01 — 10 meV [4l|, |42| enable to vary 
approximately from 0.5° to 179°, while the laser phase p is easily controlled in the range [0, 2tt). The dependence of 
\Eq (9, p)) with respect to fl and T e is analyzed considering three particular cases: 

(i) f2 3> T e : In this case the protected state becomes \Eq(tt 7 tp)) — |2) provided that — > 7r. This is achieved 
increasing the laser amplitude, since for T\ 3> T2 the lifetime of the indirect exciton state goes to infinity. This result 
is in accordance with Ref. [26| . 

(ii) ft <C T e : In the opposite scenario where 9 —> the protected asymptotic state is the ground state of the system 
\Eq (0, tp)) = |0). Since the laser amplitude is weak, the direct and indirect exciton states are not populated. 

(iii) Q = T e : For 9 = 7r/2 the protected state \E (tt/2, <p)) = (|0) — e~ lip |2)) j\[2 is a coherent superposition of 
states |0) and |2). In this last case, the control of the laser relative phase tp will enable us to implement the Deustch 
algorithm, as shown below. 

III. DEUTSCH ALGORITHM VIA DISSIPATIVE QUANTUM COMPUTATION 

The Deutsch algorithm was one of the first quantum algorithms to make explicit use of the quantum parallelism 
|15l [l7| . Such an algorithm was built to decide wheter a given binary function / : {0,1} — > {0,1} is constant 
(/ (0) = /(I)) or balanced (/ (0) ^ /(l))- Differently from the original solution to the Deutsch problem [Fa ]. 
which is probabilistic, we present here a deterministic algorithm due to Collins, Kim, and Holton |17| . Besides being 
deterministic, the approach to the Deutsch problem used in Ref. |17| is interesting because only one query to the 
oracle is made and auxiliary qubits are unnecessary. 

To implement the Deutsch algorithm we use the protected state \Eq (tt/2, p)) with the laser phase <p being or tt. 
In order to clarify the execution of the algorithm, we make the correspondence between the function domain {0, 1} 
and the states of the system {|0) , |2)} so that —> |0) and 1 — > |2). We define next the parameter e = / (|2)) — / (|0)), 
which can take the values {—1, 0, 1}. Therefore, the state \Eq (tt/2, p)) can be rewritten as 

|0)-e— |2) |Q)-(-i)*|2) |0) - (-l)/(l 2 »-/(l°» | 2 ) 
\E» (*/2, p)) = -^ = -^ = -^ . (9) 

Therefore, if the function is constant (/ (|2)) = / (|0))) then ip — 0. Otherwise, if the function is balanced (/ (|2)) ^ 
/(|0))) then p — ix. It is assumed that only the oracle has information about the function / (i). In this case the 
oracle is the laser phase p programmer. The last part of an algorithm is the readout of the solution. This is made 
here, first, by replacing the current laser field by another one with amplitude ft = Q (\/2 + l), relative phase Tp = 0, 
and frequency Zu = lu l resonant to the transition |0) — > |1). In this new configuration, \Eq (7r/2, p)) becomes the 
initial state of system and the evolved state now is ~p(t). In Fig. 1 we observe the time evolution of the populations 
Pa (t) — (i\~p(t) \i) of the states |0) (black solid line), |1) (red dashed line), and |2) (blue dotted line) considering the 

phases tp = and tp — n of the first laser. At specific times t n = nn/yfl + T£ with n = 1, 3, 5, . . ., the state of the 
system will be p (t) ~ |0) (0| if the phase of the state \Eq (tt/2, p)) is p — it and p (t) ~ |2) (2| if the phase is p = 0. 
Therefore, applying another laser pulse resonant with the transition |0) < — > |1) and observing the time-resolved 
absorption spectrum [2J] it is possible to distinguish between the states |0) and |2), provided that if the electron is 
in state |2) there will be no absorption, while if the electron is in state |0), the light will be absorbed. In summary, 
discovering wheter the phase p is or tt is equivalent to solving the Deutsch problem. 

To finish our analysis of the Deutsch algorithm, we will show that the error introduced by the fact that 
\E (ir/2,p)) is not the perfect steady state of the system is negligible. We define the quantity AP no (t n ) = 
|(0| ~p(t n ,p(t — > 00)) — p(t n , \Eq (tt/2, tp)) (E (n/2,p)\) |0)| as the difference of population in the state |0) for the 
density operator ~p(t n ) considering two different initial conditions, the approximate state \Eq (ir/2,p)} and the exact 
state pit— > 00), where the later is obtained by numerical calculation. Tables I and II show the dependence of 
APqo (t n ) on the ratio e = T2/ri for tp — and 7r, respectively. 
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FIG. 1. (Color online) Time evolution of the population of the states |0) (black solid line), |1) (red dashed line), and |2) (blue 
dotted line) considering the application of a laser field with amplitude Q. = Q U/2+ l), relative phase Tp~ = 0, and frequency 
uT = ujl- The initial state of the system is \Eq (ir/2,tp)) for (a) ip = and (b) -n. The physical parameters are Q — T e = 200 
lieV, Ti = 3 /ieV, and T 2 = 10" 4 ri. 



TABLE I. The dependence of AP o(tn) on the ratio e = r 2 /Fi for p = 0. 



APoo(ti) 

AP o(t 2 ) 

AP 00 (t3) 

APoo(t 4 ) 



e = 0.0001 



e = 0.001 



-4.9913 x 10" 
-4.9757 x 10" 
-4.9597 x 10" 
-4.9434 x 10" 



-4.9847 x 10" 
-4.9689 x 10" 
-4.9527 x 10" 
-4.9363 x 10" 



e = 0.01 



-4.9178 x 10~ 3 
-4.9010 x 10" 3 
x 10 -3 



-4.8665 x 10" 



TABLE II. The dependence of AP 00 (in) on the ratio e = F 2 /Fi for <p = n. 



APoo(ti) 

AP o(t 2 ) 
APoo(t3) 

AP 0Q (U) 



e = 0.0001 



9.9108 x 10"" 5 
9.7382 x 10~ 5 
9.5694 x 10~ 5 
9.4043 x 10~ 5 



e = 0.001 



9.8974 x 10" 4 
9.7249 x 10~ 4 
9.5563 x 10~ 4 
9.3913 x 10" 4 



e = 0.01 



9.7653 x 10" 
9.5942 x 10" 
9.4268 x 10" 
9.2361 x 10" 



IV. CONCLUSIONS AND PERSPECTIVES 

In summary, we showed the existence of a pure protected state in a system composed by two asymmetric quantum 
dots coupled by tunneling and driven by a laser field. The direct exciton state is under strong spontaneous decay. The 
asymptotic protected state is a superposition of the ground and the indirect exciton state. By controlling the ratio 
between the laser amplitude and tunneling rate it is possible to control the polar angle, while controlling the laser 
phase enables the control of the azimuthal angle of the Bloch sphere. The scheme for state protection remains true at 
low temperatures (T ~ K). For high temperatures the phonon contribution becomes important so that the phonon- 
induced dephasing will destroy the superposition of the protected state. As an application of dissipative quantum 
computation we proposed the implementation of the Deutsch algorithm in this system, which basically consists of 
distinguishing the relative phase or 7r between the states |0) and |2). Differently from the usual proposals for 
implementing quantum algorithms, which are based on unitary evolutions, here we make use of incoherent evolution 
of the Markovian reservoir. This approach is interesting because it is naturally immune to external perturbations 
and decoherence processes. A similar scheme might protect the state of two or more qubits which can be used to 
implement more sophisticated algorithms, such as Deutsch- Jozsa and Grover algorithms. 



Appendix: Effective master equation 



For the regime of parameters considered here, the master equation ([2]) has an analytical solution. However, the 
method developed in this appendix can be useful to obtain the effective master equation for a more elaborate problem. 
First we rewrite Eq. ([5} in the interaction picture according to the unitary transformation ([3]) : 



dfknt (t) 
dt 



— —i [Hint, Pint] + C-int [Pint) , 



(A.l) 



where Hint is given by Eq. ((SJ) and the dissipative Liouvillian by Eq. ([3]) with p replaced by pint . Performing a change 
of basis to the Hamiltonian eigenstates, Hi nt in Eq. (jA.lj) can be expressed as 



H int =E + (\E + ){E + \-\E_){E_\), 
and the dissipative Liouvillians C int (pint) = A (Pint) + £2 (pint) as 

£■1 (pint) = 



(A.2) 



£1 

4 



2 (9/2) [V2 cot (0/2) (\E ) (E + \ - \E ) (E_\) + \E + ) (E+\ - \E_) (E_\ 

V2cot (9/2) (\E + ) (E \ - \E_) (E \) + \E+) (E+\ - \E_) (E_\ - \E_) (E+ 
{ Pm t,\E+) (E+\ + \E„) (E_\ - \E_) (E+\ - \E+) (E_\}] 



\E-)(E + 
- \E+) <£- 



\E+)(E-\ 



Pint 



(A.3) 



and 



£2 (Pint) — —r 



sin 2 (0/2)cos 2 (6/2) \V2 cot (6/2) (\E ) (E_\ - \E_) (E \ - \E+) (E \ + \E ) (E+\) 
E+) (E+\ + \E_) (E_\-2 \E ) (E \ + \E_) (E+\ + \E+) (E_\]p int 
-V2cot (6/2) (\E a ) (E_\ - \E_) (E \ - \E + ) (E \ + \E ) (E + \) 
E+) (E + \ + \E_) (E_\-2 \E ) (E \ + \E_) (E+\ + \E+) (E_\] 



\E_)(E_\ + \E_)(E + \ + \E + )(E- 
\E_) (E \ + \E+) (E \ + \E Q ) (E+\) + 2 tan 2 



cos 2 (6/2){p mt ,\E+)(E^ 
-V^ tan (5/2) (\E ) (E_\ + \E_) (E \ + \E+) (E \ + \E Q ) (E+\) + 2 tan 2 (6/2) \E ) (E Q \}~\ , (A.4) 



where {a, b} = ab + ba states for the anticommutator. The action of a unitary transformation U (t) = exp (—iHi nt t) 
on Eq. (|A.1|) is able to remove the unitary part of its dynamics. Such procedure is interesting because the operators 
of the form \Ei) (Ej\ with i, j = {+, — , 0} for i 7^ j in Eq. (|A.1|) will oscillate quickly as shown here: 



iE+t 



UHt)\E + )(E Q \U(t) = \E + )(E \ 
UHt)\E_)(E a \U(t) = \E_){E \e^ E + t , 
UHt)\E + ) (E_\U (t) = \E + ) (E-le 2 ^. 



(A.5a) 
(A.5b) 
(A.5c) 



Since 1^ ^C Ti <C Q,T e it is possible to perform the rotating wave approximation, leading to an effective master 
equation 



dPeff (*) 
dt 



£eff (Peff) 



The dissipative Liouvillian £ e ff (.Peff) can be written in the Lindblad form as 



(A, 



C eff (Peff) = Yl Yl ~Y~ ( 20 <x,iPeff°i,i ~ °U°^Peff ~ PeffO^O^ij 



a— 1 i— 1 



(A.7) 



where 



1^ = ^ sin 2 (0/2) 




Oi,i = 


|B+){£+|-|£_)(£_| 


Ti : 2 = Tii 




Oi, 2 = 


|£-)<£+l 


Ti,3 = i\i 




Ol, 3 = 


|£ + )(£-| 


r 1)4 = ^ cos 2 (0/2) 




Oi, 4 = 


IW^-I 


Tl,5 = Ti4 




01,5 = 


\Eo)(E+\ 


r 2 ,i = £ sin 2 (9/2) cos 2 


(9/2) 


Oa,l = 


\E + ){E + \-2\E ){E \ + \E_)(E_\ 


1^2,2 = 1^2,1 




2 ,2 = 


\E-)(E + \ 


T2,3 = ^2,1 




2 ,3 = 


\E+)(E-\ 


r 2j 4 = 2r 2j i 




2 ,4 = 


cot(e/2)\E )(E + \-t&n(e/2)\E_)(E 


^2,5 = 2r2,i 




2 ,5 = 


cot (9/2) \E ) {E- 1 - tan (9/2) \E+) (E 



In Eq. (|A.7[) . the index a refers to the direct (Ti) and indirect (r2) exciton decay rates and index i enumerates the 
operators. From this effective master equation it is easy to see that \Eq) is the only eigenstate with null eigenvalue 
of operators 0\%. Although the operators 0^,% presented in C e ff (Peff) rotate the state \Eo)> the decay rates I?2,i are 
much lower than Tii. In summary, the dissipative dynamics is dictated by the direct exciton decay, ensuring that 
the only protected pure state is l-Eo). To prove the uniqueness of the protected state, we solve Eq. (|A.1|) using the 
software Wolfram Mathematica 7 for an arbitrary initial state with density matrix elements p^ (0) with i,j = 0, 1, 2 

2 

so that ^2 Pa (0) = 1. We get the expression for the density matrix elements in the asymptotic time (t — > oo): 



Poo(oo 
P\\ (oo 

P22 (OO 
P02 (OO 

Poi (oo 



r 4 (t 2 + n 2 ) r t + (r e 6 + n 6 ) r 2 

(t 2 + n 2 ) [t 2 (t 2 + n 2 ) r a + (t 4 + 2f> 4 ) r 2 ] ' 

^T2 

T 2 (T 2 + il 2 )T 1 + (T e 4 + 2tt 4 ) r 2 ' 

T 2 n 2 (r 2 + n 2 ) (Ti + r 2 ) 
(T 2 + n 2 ) [t 2 (t 2 + n 2 ) ri + (t 4 + 2m) r 2 ] ' 
r 2 (t 2 + r> 2 ) ri + (t 4 - n 4 ) r 2 



= -eTTpfi- 



(t 2 + o 2 ) [t 2 (t 2 + o 2 ) ri + (t 4 + 2r> 4 ) r 2 ] : 
P12 (00) = 0. 



Note that the elements pij (oo) are independent of pij (0). Considering the range of values for the parameters Ti, T%, 
T e , and fl used in the paper, Ti 3> 1^ and T e > fl, the matrix element pu (oo) <C poo (oo) , P02 (oo) , P22 (00). Then, 
the protected state is approximately \E (9, (p)) independently of the initial state of the system. Another way to prove 
the uniqueness of the protected state \E (9, (p)) is using theorem 2 of Ref. Q- Defining the jump operator c = |0) (1| 

and the "subspace" S = {7 \E- (9, ip)) + 5 \E+ (9, (p)) with 7, 5 € C and | 7 | 2 + \5\ 2 = l}, we observe that S ^ cS 
for every 7, 5, 9, and (p. 

The time scale necessary to the system state to be stationary is t ss — 1/ minTi^ (i = 1, . . . , 5), which depends on 
Ti, T e , and Q,. Returning to the Schrodinger picture, we observe that the unitary transformation U (t) does not affect 
the protected state, while U (t) introduces a relative phase between the states |0) and |2). 
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